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(N ! Abstract 
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cn ' We study SU(2) Lattice Gauge Theory with dynamical fermions at non- 

^ . zero chemical potential /x. The symmetries special to SU(2) for staggered 

C^ ! fermions on the lattice are discussed explicitly and their relevance to spec- 

^^ I troscopy and condensates at non-zero chemical potential are considered. Us- 
ing the molecular dynamics algorithm on small lattices we find qualitative 

J2 ' changes in the theory's spectroscopy at small and large values of fi. This pre- 
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liminary study should lay the groundwork for future large scale simulations. 

I. INTRODUCTION 

Simulations of gauge theories with dynamical fermions using the SU(2) colour group at 
nonzero density where the action is still real can be performed using standard algorithms, 
such as hybrid molecular dynamics [T[] or hybrid Monte Carlo [Q. These algorithms fail 
dramatically for QCD which has a complex action once the quark chemical potential /i 7^ 
P]. Attempts to simulate field theories with complex actions, such as QCD with a nonzero 
chemical potential, include 1. Langevin algorithms, 2. the Glasgow algorithm pioneered by 
I. Barbour 0], and 3. analytic extensions of simulations at imaginary chemical potential. 
None of these methods has succeeded. The complex Langevin algorithms suffer from uncon- 
trollable instabilities, the Glasgow algorithm requires statistics that grow at least as fast as 
the exponential of the size of the system, and simulations at imaginary chemical potential 
provide no discernable hints of transitions at real values of /x. These failures are particularly 
painful because, unlike other phenomena which are not qualitatively sensitive to the gauge 
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group, the physics of field theories at /i 7^ is crucially dependent on the gauge group. 
These points will be made clear below where SU(2) will be considered in detail. Very briefiy, 
baryons and mesons are essentially identical in the SU(2) theory, so all the special features 
of QCD at nonzero chemical potential which depend on qualitatively different meson and 
baryon spectroscopy are lost. Of even more significance to the scenarios to be discussed 
below is the fact that diquarks in the SU(2) theory can be gauge singlets while they cannot 
be singlets in SU(3). So, if we search for and find strong diquark correlations and/or con- 
densates in the SU(2) theory, these results may have no bearing on such phenomena in the 
SU(3) theory. The energetics of condensate formation is SU(3) is more challenging, both 
conceptually and computationally, than in SU(2). 

All of these points have conspired to lead researchers, both using lattice methods or 
crude analytic methods such as instanton models or mean field approaches based on gap 
equations of simplified field theories, to reconsider SU(2). This will be done below, based 
on the hope that a controlled and detailed understanding of the SU(2) theory at /x 7^ will 
give some insight into a subset of the phenomena expected in QCD. A quantitative study 
of QCD at nonzero /i must await further breakthroughs. 

The simplest scenario for QCD with /i 7^ 0, which we shall refer to as the standard 
scenario, predicts that as fi is increased, observables such as the chiral condensate (qq), which 
is non-zero in the vacuum, and quark number density n remain constant until fio ~ ^^ib/Nc, 
where mif, denotes the mass of the lightest baryon in the spectrum, whereupon a non-zero 
density of quarks is induced into the ground state, corresponding to a state of nuclear 
(or hadronic) matter. At some ^c > /^o chiral symmetry should be restored in a {first 
order) transition, resulting in a state of quark matter. Above /ic there are competing effects 
governing the quarks' influence on the gluon degrees of freedom. On the one hand, we 
expect enhanced screening due to quasi-massless virtual qq pairs as the chiral symmetries 
are restored; on the other, we expect the contribution to screening from qq pairs lying deep 
within the Fermi surface to be suppressed for kinematic reasons. 

In Figure |I] we present a tentative phase diagram for the standard scenario at zero 
temperature in SU(2) in the (m, /x) plane. The dotted line corresponds to fio{^, separating 
a vacuum phase in which n = from a phase of "nuclear matter" with n > 0. The curvature 
of this line reflects the expected behaviour mib oc ^/m, (note that the lightest baryon in SU(2) 
gauge theory is degenerate with the pion, and hence becomes massless in the chiral limit by 
standard PCAC arguments). The solid line corresponds to /ic(?Ti), which again is expected 
to scale as ^/m^. since the chiral condensate (qq) is never vanishing for m > 0, it is not 
clear a priori whether this line demarcates a genuine phase separation, or whether there 
may in principle be an end-point; this uncertainty is reflected by the question mark (the 
issue of whether there is a phase separation between "conflned" and "deconflned" phases is 
problematic with dynamical fermions). Note that in the chiral limit both critical values of 
the chemical potential tend to zero. 

Further motivation for investigation of the model comes from recent speculation that 
the ground state of dense baryonic matter may be more exotic |^ |^. On the assumption 
that at high density and zero temperature the quarks have a large Fermi surface, then it 
is natural to suppose that Bose condensation of quark pairs in the vicinity of the surface 
takes place, leading to the breaking of the global symmetry of baryon number, and hence 
superfluidity, and/or even the local gauge symmetry, leading to "color superconductivity". 



This phenomoneon would be signalled by a non- vanishing diquark condensate (qq). As 
pointed out above, for iV^ = 2 the situation is qualitatively simpler than what might happen 
in SU{3): since a SU(2) diquark condensate does not break gauge invariance, its formation 
would not lead to color superconductivity. 

The first lattice analysis of two colors QCD with non- vanishing chemical potential which 
pointed out the relevance of diquark condensates was a quenched simulation performed in 
the strong coupling limit g'^ = oo; the results were then compared with analytical predictions 
obtained from al/d expansion combined with a mean field analysis, and further confirmed by 
(nearly exact) monomer-dimer simulations [Q. For the SU(2) theory at g^ = cxo, T = 0, n ^ 
it was concluded that there is no chiral symmetry restoring transition. The simulations 
suggested that the chiral condensate (qq) actually vanishes for all /x 7^ in the limit m — > 0, 
while (qq) 7^ for /z > mT,/2, thus spontaneously breaking both the U(l)y global symmetry 
associated with conservation of baryon number and the chiral symmetry. The simulation 
results were explained by noticing that the chiral condensate rotates to a baryonic diquark 
condensate. A baryonic source term was introduced in a mean field analysis and in the 
dimer calculations but not in exact (hybrid) simulations. The spontaneous breakdown of 
the U(l)y symmetry is believed to be a direct consequence of the introduction of a chemical 
potential. A quenched study of the chiral condensate for two-colors QCD has been published 



A detailed analytic study [T^ confirms that in the instanton vacuum, at T=0, ^ > ^71^/2 
the stable phase of the two colour model is indeed a phase with diquark condensation, in 
agreement with the early lattice studies and the more qualitative discussions reviewed above. 

In section |I| we shall give a detailed outline of the global symmetries of staggered lattice 
fermions in the fundamental representation of SU(2), and show that the gauge interaction 
between quark pairs is identical to that between quark and anti-quark. In the chiral and 
zero-density limits a rotation between ground states with (qq) 7^ and (qq) 7^ can thus 
be performed at zero energy cost, and the two condensates are physically indistinguishable. 
For m 7^ 0, /i 7^ 0, the analysis must be refined; now the diquark condensate acts as 
an order parameter signalling the spontaneous breakdown of baryon number conservation. 
On the assumption that a diquark condensate forms at sufficiently high density we find it 
plausible that a genuine phase separation occurs, between a normal low density phase with 
(qq) = and a high density superfiuid phase with (qq) 7^ (as we shall see below, the (qq) 
condensate which seems most likely to form is gauge invariant, and does not give rise to 
superconductivity in this case). This is shown as a solid line in the tentative phase diagram 
of Figure |^; the curvature of the line this time reflects the fact that at large Fermi energies 
asymptotic freedom should weaken the interaction between quark pairs and thus favour the 
chiral condensate as m is increased. Once again the onset for thermodynamics is expected 
to occur for fio ex m^, and this is shown by the dotted line. 

By further increasing 11 > ^^ — it^-k/'^ diquarks could begin to populate the vacuum 
producing the ordinary cold nuclear matter phase of SU(2) gauge theory. If the dynamics 
favor pairing, they could condense. Of course, this phenomenon would have to compete with 
the free quark phase. Computer simulations which are sensitive to relatively small energy 
differences could, in principle, find the separation between ^o — 'm.,^/2 and /ic and assess the 
physical phase between the solid and dotted line in the phase diagram outlined in Figure |^. 

In this paper we present an analysis of the symmetries of the staggered lattice fermion 



action. We identify Goldstone modes appropriate to the staggered fermion lattice action 
in the regimes of low and high baryon density, which enables a tentative classification of 
the spectrum of light excitations in each regime. We then present numerical results for 
unquenched simulations of two colours QCD at /? = 1.3 on small lattices, which suggests 
that for Nc = 2 the rotation of the chiral condensate into a diquark condensate observed at 
infinite coupling is likely to occur also at this coupling, and we establish the level ordering 
in the particle spectrum in the high energy phase. 

II. ACTION AND SYMMETRIES 



We will outline plausible patterns of symmetry breaking for SU(2) lattice gauge theory 
with fermions in the fundamental representation of the gauge group. We emphasize the 
distinction between the model with continuum fermions in a pseudoreal representation (such 
as the 2 of SU(2)) which has the symmetry breaking pattern SU(2iVj) -^ Sp{2Nf) in contrast 
to the lattice model with staggered fermions which, as we shall outline below, has the 
unorthodox pattern SU(2A^/) — > 0{2Nf). We also consider the consequences of introducing 
a chemical potential corresponding to non-zero baryon density. 

Let us start with the kinetic term in the lattice action for a gauged isospinor doublet (2) 
of staggered fermions x? X- For clarity we will consider A^ = 1 flavors, although this will be 
generalised when the numerical simulations are discussed. 
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The [/'s are 2x2 complex matrices acting on color indices. The chemical potential enters 
via the timelike links in the standard way. It is straightforward to rearrange ( p.l|) , using 
the Grassmann nature of X)Xi the fact that ti^{x) = (— 1)^o+"=^a'-i = ri^{x ± /t), and the 
property T2U^T2 = U*, where T2 is a Pauli matrix: 
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(2.2) 
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The sum in ( |2.2| ) is only over even sites, ie. those defined by e{x) = (^—iyi+^2+x3+x4 _ _j_]^^ 
and the subscripts emphasise that the X fields are only defined on even sites, and the X on 
odd. 



Now, in the original action ( |2.1|) there are two manifest global U(l) symmetries: 

X ^ e'"x X ^ Xe-*" : X ^ e'"'x X ^ Xe*"^ (2.4) 

the first (U(l)y) corresponding to conservation of baryon number and the second, which 
holds only in the chiral limit m = and which we will call U(l)^, corresponds to a form of 
axial charge conservation. In ( |2.2| ) we can see that in the chiral limit at zero density (ie. for 
fi = m = 0) both are subsumed in a larger U(2) symmetry (because exact global symme- 
tries are non-anomalous in lattice formulations, the full symmetry is U(2); in a continuum 
approach the axial anomaly would reduce it to SU(2)): 

Xo ^ VXo Xe ^ XeV^ V G U(2). (2.5) 

By introducing a non-zero chemical potential we reduce the U(2) symmetry of the lattice 
action to U(l)y ® U(l)^, and by introducing a bare mass we reduce it to 11(1)^. This will 
be important when we identify and compare the Goldstone modes present in the low and 
high density phases. 

A. Qualitative Discussion of Symmetries and Spectrum 

1. ^ = 0.0 

As anticipated in the Introduction, when /i = quarks and antiquarks belong to equiva- 
lent representations: there is a symmetry transformation which turns a chiral condensate into 
a diquark one. In other words, there is only one discernable condensate, < qq >'^ + < qq >^, 
whose orientation in the chiral sphere is selected by the explicit breaking term. 

Recall that the extra quark-antiquark symmetry enlarges the continuous global symmetry 
of the lattice Action from U{f) x U{f) to U{2f), while the subgroup which leaves the mass 
term < qq>^ + < qq >^ invariant is 0(2f) [llj. By applying a generic 0(2f) transformation 
to the mass term XX we identify the following set of operators which might form condensates: 

scalar xiXi + X2X2 pseudoscalar ^(xiXi + X2X2) (2.6) 

scalar diquark X1X2 — X2X1 pseudoscalar diquark £^(xiX2 — X2X1) (2.7) 

scalar antidiquark X1X2 — X2X1 pseudoscalar antidiquark £(xiX2 — X2X1) (2.8) 

where the lower index labels colour. The first line displays the usual pseudoscalar and scalar 
condensates. The second (third) line corresponds to diquark (antiquark) condensates, scalar 
and pseudoscalar. This simple minded quantum number assignment is done by considering 
that quark - quark and quark-antiquark pairs have opposite relative parity, and will be 



confirmed by the rigorous analysis below. See also the interesting discussions in |T^ . 

Consider now quark propagation from a source at to the point x. The propagator Gij 
{i,j color index) is an SU(2) matrix: 

G,= (_; M (2.9) 



structed from Gij as foUows:- 




pion 


tiGG^ = {a^ + b^) 


scalar meson 


etiGG^ = eia" + b^ 


scalar qq 


detG = (a^ + b^) 


scalar qq 


detGt = (a^ + b^) 


pseudoscalar qq 


edetG = eia" + b^ 


pseudoscalar qq 


edetG^ = e{a^ + b^ 



By taking its correlations in the appropriate sector of quantum numbers we form meson 
and diquark operators. We shall limit ourselves to the local sector of the spectrum and 
focus on the zero momentum connected propagators of the scalar and psudoscalar mesons 
and diquarks. The scalar meson propagator will thus be an isovector, which we will call S, 
following QCD notation. 

The meson (qq) and diquark [qq) and antidiquark {qq) propagators at /i = are con- 

(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 

The notable feature of the propagators at yU = is the exact degeneracies of the pion, 
scalar qq and scalar qq and of the scalar meson, pseudoscalar qq and pseudoscalar qq. This 
can also be understood as being due to quarks and anti-quarks having opposite intrinsic 
parities, whereas the qq and qq interactions due to gluon exchange are identical. We then 
identify two orthogonal directions in the chiral space: 
a) TT -scalar diquark - scalar antidiquark 
h) S - pseudoscalar diquark - pseudoscalar antidiquark 

2. fi^O 

When // 7^ the symmetry is reduced to that of staggered fermions, i.e. ^a(I) x Uv{l) 
P] , however the possibility of diquark condensation Q could lead to unusual patterns of 
chiral symmetry. 

The condensate will tend to rotate in chiral space as /i increases, rotating into a purely 
diquark direction for large /i. However, as /i increases and the symmetry in chiral space 
is reduced, as will be discussed below, the rotation will no longer be "trivial" (as it would 
be at zero chemical potential) and the new vacuum would be physically distinct from the 



original one. See also the interesting comments in |T0| . 

Consider the Dirac kinetic operator p for staggered fermions. At /i = the relation 
p'^ = —p holds. For /j, ^ the fact that e^ {g~^) multiplies the forward (backward) 
timelike gauge links in the lattice action means that p"^ ^ ~Pi because in the matrix p'^ , 
g-/x ^g/x^ j^Q^ multiplies the forward (backward) links. Consequently the propagators G and 
G^ are no longer trivially related as they are at /i = 0. The composite propagators can be 
computed as follows |T^ : 

pion tiG{fi)G\—n) 
scalar qq detG(/i) (2.16) 



^ The conditions leading to the Vafa Witten theorem do not hold here. 



with analogous formulas for the scalar meson and pseudoscalar qq. 

Note for /z 7^ 0, that although the diquark Goldstone modes in ( 2.24 ) contain both xx 
and XX operators, we expect the physical states in the high density ground state to be 
formed predominantly of quarks rather than anti-quarks. We see from ( p.l6|) that for /i 7^ 
the pion is no longer degenerate with the scalar qq; correspondingly the degeneracy of the 
scalar meson with the pseudoscalar qq is also broken. 

Finally, consider the "baryonic" pion |TB[ which is believed to be responsible for the 
failure of quenched 3 colours QCD (or indeed for any treatment with a real path integral 
measure proportional to detM^M). The propagator would be: 



baryonic pion trG'(/i)G'^(/i) 



(2.17) 



We see that for the 2 colours theory a consequence of the fact that the quark propagators 
are SU(2) matrices even for /x 7^ is that the baryonic pion is degenerate with the scalar 
diquark via eqns. ( p. 101 , ^1121) : in SU(2) at nonzero chemical potential it is still possible to 
express a diquark propagator as a positive definite quantity. 



B. Rigorous discussion of symmetries and spectrum 



At zero density (/i = 0), we expect the model to display spontaneous chiral symmetry 
breaking, with a chiral condensate of the same form as the mass term mJ2xX{^)x{^)- I^i 
terms of the X,X fields, introduced in (E.3|), 



/C/C r) 



X. 



1 



r2XT + K 



tr 



1 



rsX, 



(2.18) 



Of course, the condensate breaks the global symmetry ( |2.5| ): the residual symmetry is 
generated by the subgroup which leaves invariant the symmetric form ( J J ) , whose most 



general element is 



which generates U(l)y. We thus identify the pattern of 



chiral symmetry breaking as U(2)^U(1) for /i = 0. 

For /i = the number of broken generators is 3 (since dim U(2)=4), so we expect the 
spectrum of the model to contain three massless Goldstone modes. These can be found by 
considering infinitesimal rotations of the condensate ( |2.18 ) by V^ = 1 + i5\, with A one of 
the U(2) generators {1, Tj}, and identifying the mode as the coefficient of 0{5). The results 
are 
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(2.19) 



As expected, the first of these modes is the familiar pion x^Xi which is a pseudoscalar since 
it is odd under the lattice definition of parity: 



X = {xo,Xi,X2,X3) 



X[X) 



4)<+4^(^/ 



X 



X[x, 



(xo,l -Xi,l -X2,l -X3) 
^ (-1)<+-3;^(X'). 



(2.20) 



The second two, however, are gauge invariant hnear combinations of scalar diquark xx ^^^ 
XX states. The 1 indicates that rotations generated by T3 leave the condensate invariant. 
These states carry baryon number and although they would be true Goldstone modes only 
in the chiral limit and at /x = 0, by continuity we expect them to be light states in the low 
density regime with masses vanishing as y/m. For /x 7^ of the three states identified only 
the pion remains a Goldstone mode. The pion is generated by 1 and is associated with 
spontaneous breaking of ^7(1)^. The U(l)y rotation generated by t^, as discussed above, 
leaves the condensate invariant. 

The symmetry breaking described here is notable because it is distinct from the 
SU(2A^j) — i> Sp(2Nf) predicted for continuum fermions in a pseudoreal gauge represen- 
tation such as the 2 of SU(2) |T^. This discrepancy between lattice and continuum has 
been known for some time [|ri|] ||T6[. For a complete classifications of the possible patterns 
of chiral symmetry breaking see |T^ . For patterns of chiral symmetry breaking the roles of 



fundamental and adjoint representations appear to be reversed on the lattice |TB|. 

At high density (/i ^ m^/2), let us postulate that a large Fermi surface will promote the 
formation of a diquark condensate. Which condensate forms is a dynamical question, since 
in principle many diquark wavefunctions can be written down. It is instructive to consider 
a "maximally attractive channel" (MAC) for the formation of the condensate. We proceed 
by making the following assumptions about the MAC condensate: 

(a) The condensate is gauge invariant. 



(b) The condensate is invariant under lattice parity ( |2.2CI| ) and exhibits no preferred spatial 
orientation. 

(c) The MAC condensate is as local as possible in the x? X fields. 

(d) The condensate wavefunction is anti-symmetric under exchange of fields. 

Of these four assumptions, only (d) should be considered inviolate |^ [H. For isospinor 
fermions, however, it is straightforward to write down a local condensate qq2 which satisfies 
(a) - (d): 



9?2 



1 r 



X ''ix)T2Xix) + xix)T2X ''(x) 



(2.21) 



The relative + sign between the two terms is not arbitrary; this can be seen by considering a 
two-point function {xx{x)xx{y)) , which in a vacuum with a diquark condensate will cluster 
at large separations into {xx){xx) IHI H3- We know that the scalar diquark propagator is 
proportional to the determinant of an SU(2) matrix (see eqs.( p.l^ - p.l5| ) above), and hence 
positive definite; therefore the xX ^^^ XX condensates must form with the same sign. 
In terms of the X, X fields the condensate is written 



qq2 = ^ 



X, 



^ _Jr2Xr + Xr^^ 



T-2X„ 



f2.22) 



The fascinating feature of this SU(2) model is that the diquark condensate ( p.22|) can be 
obtained from the chiral condensate ( p.l^ ) by an exphcit global U(2) rotation: 

Therefore the pattern of symmetry breaking for m = fi = remains U(2) —>■ U(l), but the 
unbroken symmetry is no longer U(l)y. The chiral condensate is equivalent to the diquark 
condensate due to the symmetry transformation relating them. There is only one discernable 
condensate whose orientation in the chiral sphere is selected by the explicit breaking term. 
In lattice simulations we select the chiral condensate (qq) because the bare mass m ^ 0. In 
principle by setting the bare mass to zero and including a diquark source term we would 
obtain a diquark condensate (qq) at fi = 0. 

For non-zero density one might naively expect that as /i increases with respect to m 
the condensate gradually rotates from a chiral {qq) to a diquark {qq). Actually the true 
behaviour is more subtle, because once m ^ 0, fi y^ the rotation is no longer trivial; the 
vacuum at large n must be distinct from that at /i = 0. This fact is reflected by different 
numbers of Goldstone modes in high and low density regimes; we shall argue that there is a 
phase separation between low and high density, with the diquark condensate acting as order 
parameter. 

We can check this by analysing the Goldstone modes obtained by U(2) rotations of the 
diquark condensate (|2.2 



(fT2(): 






1 => x^''T2ex + XT-2£X*'' 




Ti ^ 1 




^2 ^XX 




Tz => X^^'riX - Xrax*'' 



(2.24) 

Just as before, in the chiral limit for /x = 0, there are three Goldstone modes, a scalar qq 
state, a scalar qq state (which is orthogonal to the condensate), and a pseudoscalar qq state. 
The spectrum of massless states has identical quantum numbers to those of ( |2.19| ) in this 
limit. The pseudoscalar gg, i.e. the pion, is no longer a Goldstone particle in the high 
density vacuum. For yU 7^ we retain two of these Goldstone modes: firstly the pseudoscalar 
diquark generated by 1 and associated with spontaneous breaking of U(l)^ which is a 
pseudo-Goldstone for m 7^ 0; secondly the scalar diquark generated by t-^ and associated 
with spontaneous breaking of U(l)y. This second state remains an exact Goldstone mode 
even once m 7^ 0, since there are no terms in ( p.l|) which explicitly break this symmetry. 
Therefore we expect a level ordering (light — ^ heavy) in the high density regime of scalar 
diquark (massless), pseudoscalar diquark, scalar meson, pion. We will see below in sections 
[rV| and |V| to what extent these predictions are borne out in the simulations. 

Recall that away from the chiral limit, in the low density regime but for /x 7^ we have 
no exact Goldstones, but only three light states, the pion and two scalar diquarks. The fact 
that the low and high density regimes are characterised by different numbers of massless 
modes indicates that they should be separated by a true phase transition at some ^ = fie- 
This is indicated in the proposed phase diagram of figure ^ The Goldstone count in the 
diagram is as follows: zero in the region to the right of the solid line; one in the region to 

9 



the left of this hne where {qq) ^ 0; two in the chiral hinit along the /x-axis; finally three at 
the U(2) symmetric point at the origin. 

Finally, we should discuss the effects of using the hybrid molecular dynamics [|I| (or 
hybrid Monte Carlo 0) algorithm to simulate the model. The functional integral measure 
contains a factor det M^ M^ with M the fermion matrix for one flavor of staggered fermion. 
Since detM''' = detM* = detr2Mr2 = detM, then we know the simulation describes 2N 
identical flavors of staggered fermion. It can be shown that detM is positive definite for 
/i 7^ |21[] therefore one could in principle "take the square root" by setting A^ = | in a 



hybrid algorithm. To reduce the number of flavors in a simulation at zero chemical potential 
using staggered fermions, standard practice is to use even-odd partitioning. This relies on the 
equality of the fermion determinant on the even sublattice with that on the odd sublattice 
which means that we require det M^ = det M, a relation which no longer holds for /i 7^ 0. 
Thus we are constrained to a minimum of Nj = 8 continuum flavours in our simulation 
(recall that the number of physical flavors Nf = AN for staggered lattice fermions). This 



should still result in an asymptotically-free theory, since Nf < -^Nc, however, the two- 
loop coefficient of the beta-function is now positive since Nf > 34N^ / {13N^ — 3), which 
means that the model should exhibit non-trivial fixed-point behaviour in the far infra-red 
P^ . The question of confinement /deconfinement is thus further complicated. A recent non- 



perturbative calculation by Appelquist and Sannino |Q suggests indeed that the onset of 
the conformal phase is at Nf > 3.9Nc. That would place our 8 fiavours 2 colours results 
in the conformal phase; however, on the lattice volumes considered here we have seen no 
indication of an 'exotic' behaviour. 

The symmetries of the lattice model for N = 2 fiavours of staggered fermion will be 
very similar to the N = 1 case discussed above because (i) the mass term and hence the 
chiral condensate will be fiavor symmetric, and (ii) the diquark channel we have analysed 
on the assumption that it is MAC must also be fiavor symmetric due to the Pauli exclusion 
principle. Therefore the whole analysis should go through virtually unchanged, with the 
resulting pattern of symmetry breaking for /i = of U(4) —>■ 0(4), with 10 associated 
Goldstone modes. Of course, all our measurement algorithms to date work in the single fiavor 
sector, therefore we should only "observe" three Goldstone modes. So the special features 
of the SU(2) gauge theory arise from the fact that the quarks and the anti-quarks belong 
to equivalent representations which enlarges the continuous global symmetry of the lattice 
action describing N fiavors of staggered fermion from U(iV) x U(iV) to U(2iV). The subgroup 
which leaves the mass term invariant is 0(2A^). In general there are 4iV^ — N{2N — 1) = 
2N'^+N broken generators associated with N"^ Goldstone mesons, and N'^+N = 2N{N+l)/2 
Goldstone baryons and antibaryons. It is instructive to consider what will happen as the 
flavor symmetries enlarge in the continuum limit, when we expect the model to describe 
Nf = 4N quarks. It is currently unclear whether the lattice model in this limit will persist 
in the non-standard SU(2iVj-) — >■ 0{2Nf) (note the axial anomaly should reemerge in this 
limit), or assume the orthodox continuum pattern SU(2iVj) — > Sp{2Nf) [|15[] . 

III. SIMULATIONS 

In this initial study of two colours QCD at non-zero chemical potential we used a standard 
hybrid molecular dynamics algorithm P|,P^ . Since we cannot use even-odd partitioning when 
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/i 7^ we were constrained to the minimal value of 8 continuum flavours. 

We simulated at /x 7^ at couplings: /? = 1.3 where the system was in the phase of 
broken chiral symmetry at /i = 0. Most of the simulations were on 6^ lattices however we 
studied the mass spectrum on a 6^ x 12 lattice at /3 = 1.3 with m = 0.07 for /i = 0.0, 0.2, 0.4. 
Note that although the fermion determinant of the SU(2) theory is positive deflnite for 
yU 7^ 0, it is likely that value of the fermion determinant has large fluctuations in magnitude 
when the chemical potential is switched on. In support of this conjecture we found that the 
simulations numerically intensive for /i > 171^^/2. 

For chemical potential smaller than m^/2 simulations were not problematic and in this 
case we chose dt = .02 and ran for 5000 sweeps, after discarding 1000 sweeps for thermal- 
ization. 

For /i > m7r/2, we observed several instabilities. Typically, the system was crashed after 
a few tens of iterations. To combat this we refreshed the system frequently (every 5 sweeps 
for the gauge flelds and 3 for the fermions). In this way we could run for thousands of 
iterations with dt = .02. With the exception of /i = .8 , dt = .05 worked also. As flnal 
strategy, we thermalized with dt=.05 for 500 iterations, subsequently reduced dt to .02 , 
reequilibrated for 300 sweeps proceeded to measure for 5000 sweeps. 

To monitor any possible bias introduced by the frequent refreshment we also continued 
our runs at /i = .6 with a pure molecular dynamics evolution, which revealed no problem 
for a few thousands steps (therefore we stopped it). 

Clearly, we cannot exclude ergodicity problems. 

In Figs. ^ we show the chiral condensate as a function of bare mass m for a range of 
values of chemical potential. The system is clearly in the broken phase at /i = 0.0 and the 
chiral condensate decreases quite sharply as /i is increased. 

In the standard scenario, we expect that as /x is increased beyond m^/2 the ground state 
populated by quarks will be favoured over the vacuum state; therefore chiral symmetry may 
be restored for /i > ■m^/2. We found that at /3 = 1.5 for ?n = 0.1, the pion mass m.,^ ~ 0.8 
while at /? = 1.3 with m = 0.07 we flnd m.,^ ~ 0.6 

Behaviour consistent with the standard scenario is reflected in the Figs. ^ where the 
chiral condensate is plotted as a function of m. At (3 = 1.3 and at // = 0.2 there is a strong 
m-dependence in the value of the chiral condensate. The chiral transition seems signiflcantly 
sharper at m = 0.07 than it is at m = 0.05 as we would expect. See the spectrum section [V| 
for a discussion of these results. 



IV. PARTICLE SUSCEPTIBILITIES 

Measurement of the meson and diquark susceptibilities will allow us to assess the effects 
of chiral symmetry restoration and, to some extent, explore the level ordering in the particle 
spectrum [^] which should reflect properties noted in the earlier discussion of symmetries 
where we identifled Goldstone modes in high and low density phases. 

In terms of the timeslice propagator P(t) = J2xG (x, t), the susceptibility, X) is deflned 
by 



X = j:Pit). (4.1) 
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On the assumption that a particular channel is dominated by a single pole, then x can be 
written Z/M'^^, where M^g is a bound state mass and Z is a wavefunction renormalisation 
constant. If we make a second assumption that the Z's are relatively insensitive to the 
quantum number of the state, then susceptibility measurements yield information about 
level ordering without the need for any propagator fitting, and may be a more reliable guide 
on the small lattice volumes presented here (see however the comments below). 

We have measured susceptibilities of the scalar mesons and diquarks, taking into account 
contributions from connected diagrams only, and plotted the quantity mx which we expect 
(from the Ward identity) to be equal to the chiral condensate in the case of the pion suscep- 
tibility. Although the realization of PCAC in this model has yet to be worked out in detail, 
it is reasonable to assume that in each channel the susceptibility gives a direct measure of 
the associated condensate. 

The degeneracies of the meson and diquark sectors are broken for any /i 7^ 0. The pion 
which is a pseudo-Goldstone for // < m^/2, becomes heavy in the dense phase whereas the 
pseudoscalar qq state which is heavy in the low density phase, becomes light in the high 
density phase. Recall from section ^ that we expected the scalar diquark to be an exactly 
massless mode in the high density phase. Unfortunately it is impossible to infer the level 
ordering directly from the susceptibilities, since the possible occurrence of a plateau prevents 
it (see the discussion in Sec. |V]). 

The level ordering we infer (lightest to heaviest) from Fig. § for the susceptibilities is 
scalar diquark, pseudoscalar diquark, and a degenerate pion and delta susceptibilities. 

Degeneracy of the pion and 5 is usually interpreted as a signal for the restoration of ^/^(l) 
symmetry. It would be natural to infer from that a complete suppression of the topological 



fluctuations induced by the ordering effects of the chemical potential |26 



Rigorous arguments predict that, indeed, the scalar diquark propagator should always 



be the larger one ||2^. However, it is not possible to infer from this a mass level ordering 
since the scalar propagator has a disconnected component, and because of the multiplicative 
renormalization of the wavefunction. For these reasons we have performed a mass spectrum 
calculation which will show that the particle content of the scalar and pseudoscalar diquark 
propagators is indeed the same and we postpone a detailed comparison of the numerical and 
analytic results to the spectrum section. 



V. PARTICLE SPECTRUM 

To extend and verify the information obtained from the susceptibilities we have per- 
formed particle spectrum measurements on a lattice with extended temporal direction. We 
simulated at /3 = 1.3, m = 0.07 on a 6^ x 12 lattice. The qualitative trends revealed in the 
susceptibility measurements are the fact that we have a light pion and scalar diquark in the 
low density phase and two light diquark states in the high density phase. We expect the 
pion to become heavy and degenerate with the connected contribution to the scalar meson 
(called 6) for fi > mT,/2. 



Finite density spectroscopy analysis has been introduced and discussed in |]I4|. The 
following symmetries should hold true in the ensemble: 

G{t) = G{T - t) for mesons (5.1) 
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G(t, jj) = G{T — t, —jj) for diquarks (5.2) 

The particle propagators for /i = 0.0 and for /i = 0.4 are shown in Fig. 6. At /x = 
the pixon and scalar diquark are degenerate light states while the pseudoscalar diquark and 
scalar meson are degenerate but comparatively heavy. The propagators are symmetric in 
time at ;U = 0. However at /x = 0.4 the chemical potential induces an asymmetry in the 
propagators by favouring forwards over backwards propagation. As expected from the sus- 
ceptibilities, the pion is now a heavy state which is approximately degenerate with the scalar 
meson. We also note (see Figs. 6 and ^ ) that at /x = 0.4 Gscaiar qq(i) = Gpseudoscaiar qq(i) +P, 
where p is a constant. This means that the pole contents of scalar and pseudoscalar diquark 
propagators are identical, i.e. the scalar and pseudoscalar diquarks are degenerate. We 
found that p increased as /x was increased, this trend being more clear in the scalar channel 
where condensation is expected (Fig. ^. If this trend were to persist for a range of lat- 
tice volumes then it could be intepreted as evidence for condensation in the scalar diquark 
channel with p oc \{qq)\'^- A similar analysis has been applied to simulations of the 2+1 
dimensional Gross-Neveu model [[19|; stable plateaux were indeed observed, but the volume 



scaling of the constant p did not show the expected behaviour. To confirm this scenario one 
would have to implement an explicit diquark source term jqq and extrapolate the source j 
to zero: this work is in progress |2^ . 



The dashed lines in Figs. 6 are the results of the fits to 

Gra{t) = a(e(-'"*) + e'"^^-*)) + (-l)*6(e(-"i*) + e^^^^"*)) + p (5.3) 

for the mesons and 

G'fc(t) = ae"™^* + 6e'"''(^-*) + p (5.4) 

for the diquarks. 

For the mesons p should be zero and m and mi correspond to masses of states with 
different lattice parity. Since the mesons do not "feel" the effects of the chemical potential 
the masses of the forward and backward propagating meson states are identical. Although 
the inclusion of the oscillating term was mandatory for a good fit, very often only the 
fundamental amplitude and the first mass were meaningfully determined. For /i 7^ we 
expect iTif 7^ mij in the diquark channels, reflecting the different backwards and forward 
propagations in the dense medium. 

The particle masses for the mesons and for the diquark states were measured at /x = 
0.0, 0.2, 0.4 and are plotted in Fig. |^. The pion mass increases, while the scalar meson 
mass decreases as /x is increased and the two particles become approximately degenerate at 
/i ~ 0.4 (at zero density vn.„j1 ~ 0.3). The scalar diquark which is degenerate with the pion 
at yU = 0.0 is already lighter than the pion at // = 0.2 and by /x = 0.4 the pseudoscalar and 
scalar diquarks are apparently considerably lighter than the pion and delta meson. 

In conclusion the ordering is as follows: two light (see however the comment below), 
degenerate diquarks, two heavier, degenerate pion and delta 

This ordering is not quite what was predicted in section ^ and several comments are in 
order. First, recall that the predictions of section |I| are for the a meson, while the numerical 
results are for the b meson. It is certainly possible that the a is lighter. 
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Next, the operator ( |2.12|) , being manifestly real, may not project onto the true Goldstone 



mode, which corresponds to fluctuations in the phase of the condensate. It is interesting 
to note that for diquark condensation, unlike chiral condensation, there is no spacetime 
quantum number (such as parity) which serves to distinguish transverse (ie. Goldstone) 
from longitudinal (ie. Higgs) excitations. 

A speculative comment is in order. A naive approach to calculating masses in the 
presence of a chemical potential would suggest that the zero density masses would be shifted 
downwards by n x /x, where n is the number of quarks in the state under consideration. This 
shift merely reflects a change of the reference level of the energy, therefore does not take 
into account any new dynamics associated with propagation through a dense medium. To 
appreciate the true dynamical content we plotted, together with the masses themselves, the 
reference energy -^(/i) = -E(O) — 2/i. We do indeed observe that the change of the masses is 
genuine and not merely a simple shift in the energy scale. However, if we add a shift n/x to 
our spectrum results, the data suggest that all the four masses are then nearly degenerate. 



Note that preliminary calculations in a model field theory suggest precisely this result [p^j .j 
Finally, if we estimate the constituent mass from the p propagator, as seems reasonable, 
we would convince ourselves that these states are still bound as p ~ 1.8 > vr ~ 1.2. In 
conclusion, the results for the n-S, and the scalar-pseudoscalar diquarks, masses suggest 
degenerate, composite, massive, but rather hght particles. 

VI. SUMMARY AND OUTLOOK 

In this study of two colors QCD at non-zero chemical potential with dynamical fermions 
we have outlined in detail the symmetries of the lattice model and identified the Goldstone 
modes in the low and high density phases incorporating the possibility of diquark conden- 
sation for /i > mjr/2. In our preliminary simulations we have established the level ordering 
of the particle spectrum from the susceptibilities and via direct measurements for several 
values of /i. The analytic results predict Goldstone states which we have not observed in 
the simulations. The simulation results indicate four degenerate states at large /i, when the 
data are referred to the same energy scale, or two couples of degenerate particles if we do not 
change the energy reference. As discussed in the text, it is conceivable that none of the oper- 
ators we have implemented are suitable for measuring the Goldstone modes predicted by the 
analytic calculations. We have noted exact degeneracies at /i = between pion and scalar 
diquark; scalar meson and pseudoscalar diquark. The pion which is a pseudo-Goldstone in 
the low density phase becomes heavy in the dense phase whereas the pseudoscalar diquark 
is heavy in the low density phase but becomes light in the high density phase. For /i > m^/2 
the pion and scalar meson 5 are approximately degenerate. It would be very interesting to 
measure the a which would require the disconnected contributions. It is certainly possible 
that the sigma and delta masses are different. 

It is possible that the rotation of chiral condensate to baryonic condensate first noted in 
the infinite coupling regime J^ persists in the intermediate and weak coupling regimes and 
our symmetries analysis supports this. The fact that the scalar and pseudoscalar diquark 
states are found to be light for /i > m^/2 is also suggestive. In an extension of this work 
we will implement an explicit diquark source term in a hybrid Monte Carlo simulation and 
investigate the effects of extrapolating the source to zero 

14 



Finally, we notice that while the diquark observables introduced for studying conden- 
sation phenomenon display an interesting and peculiar behaviour, there is nothing in the 
behaviour of the conventional observables discussed here - - chiral condensate, meson spec- 
trum - which requires going beyond the standard scenario. The thermodynamics of this 



model is considered in a companion paper |29 
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FIG. 1. Tentative phase diagram in the {m,n) plane for the standard scenario 
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FIG. 2. Tentative phase diagram in the (m, /i) plane when a diquark condensate {qq) is con- 
sidered 
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